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A three-phase unbalanced power flow method can provide a more realistic 
scenario of how distribution networks operate. The backward/forward sweep- 
based power flow method (BF-PF) has been used for many years as an 
important computational tool to solve the power flow for unbalanced and 
radial power systems. However, some of the few available research tools 
produce many errors when they are used for network reconfiguration because 
the topology changes after multiple switch actions and the nodes are 
disorganized continually. This paper presents a modified BF-PF for three- 
phase unbalanced radial distribution networks that is capable of arranging the 
system topology when reconfiguration changes the branch connections. 
A binary search is used to determine the connections between nodes, 
allowing the algorithm to avoid those problems when reconfiguration is 
carried out, regardless of node numbers. Tests are made to verify the 


Power flow usefulness of the proposed algorithm in both the IEEE 13-node test feeder 
Radial Network and the 123-node test feeder, converging in every run where constraints are 
accomplished. This approach can be used easily for a large-scale feeder 
network reconfiguration. The full version of this modified backward/forward 

sweep algorithm is available for research at MathWorks. 
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1. INTRODUCTION 

Distribution networks are unbalanced systems in which loads are not precisely balanced through the 
network [1]-[3]. This feature rejects conventional balanced power flow methods, such as Newton—Raphson 
and Gauss-Seidel, to converge in order to find a solution for the non-linear power flow problem [4]-[6]. 

BF-PF has been used for several years to solve the radial power flow problem and is based on both 
the current and voltage Kirchhoff’s laws. This method was initially implemented for solving balanced radial 
distribution networks [7]. Nevertheless, the most commonly found distribution networks are unbalanced 
because most residential and commercial users use single or two-phase loads for home appliances and office 
equipment. On the other hand, most industrial users request three-phase feeders due to their load 
demands [8], [9]. 

Wang et al. presented an unbalanced BF-PF for real distribution systems [10], [11]. Diaz et al. 
developed a direct BF-PF for solving load power flows in AC droop-regulated microgrids [12]. Cun et al. 
performed a fuzzy approach to evaluate the power flow through the BF-PF method [13]. Kumar and Jain 
used the BF-PF algorithm for calculating both power losses and voltage regulation in radial systems with 
distributed generators and voltage regulators [14]. Nimal Madhu et al. incorporated a droop control to the DC 
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equivalent power flow method for distribution and low voltage systems [15]. Xue et al. implemented an 
unbalanced three-phase distribution system power flow with distributed generation using the affine arithmetic 
self-validation method. An active distribution power flow analysis was performed using asymmetrical hybrid 
technique in [16]. Issicaba and Coelho proposed some improved methods inspired on the classical Cespedes’ 
load flow method, improving its efficiency in time response [17], [18]. 

Other authors have used the BF-PF for calculating load flow in distribution 
systems [19], [8], [9], [20]-[26]. With those improvements, BF-PF has become a useful tool for solving 
distribution network problems. Most of the approaches found in the literature and available for research for 
the power flow problem, consider a consecutive order in the node numeration. A downstream node should 
have always a higher number than its feeder bus. As an example, node number 2 is always fed by node 
number 1, never the contrary. 

This paper presents a modified backward/forward sweep-based power flow method for running the 
power flow of unbalanced radial networks that is suitable for the reconfiguration of distribution systems. 
The changes apply a binary search process to search the nodes, regardless of their numeration and 
organization. In that way, node number n could be fed by node number m. This approach allows researchers 
to make changes in the basic topology on any three-phase unbalanced radial distribution system. Considering 
the above, a feeder network reconfiguration process can be achieved even if the node numeration is not 
identified by a strict ascendant order. 

The rest of the document is organized as follows. Section 2 presents a brief explanation concerning 
the standard BF-PF method and its implementation on IEEE 13-node test feeder. Section 3 explains 
improvements developed in this paper, considering the binary search process for identifying the topology of 
the system. Section 4 presents the modified IEEE 13-node test feeder and the IEEE 123-node test feeder to 
make a reconfiguration. The evaluations in Section 5 provide results and discussion that guarantee the 
proficiency and usefulness of the proposed method. Finally, some conclusions are presented in Section 6. 


2. STANDARD BF-PF 

The BF-PF is derived from both the voltage and current Kirchhoff’s laws [27], [28]. In Figure 1, we 
present a flowchart with the explanation of the steps of a standard BF-PF, which is commonly used for 
solving the power flow problem. We have divided this flowchart into three main sections: the input values, 
the process, and the output values of the algorithm. 

The inputs are given according to the distribution network test case and, for this test, we considered 
a total of seven inputs. The first input data correspond to the line segment, which provides information about 
the connections between the nodes of the radial distribution network. The second input is the voltage matrix 
that contains the voltage magnitudes and angles per phase at each node. The third input is a load matrix that 
contains information of the real and reactive power loads at each node. The fourth and fifth inputs are the 
identification of the main feeder as the slack node and the sources with voltage regulation as the PV nodes to 
maintain the voltage magnitude at each iteration. The sixth input is the impedance configuration matrix, 
given by the distribution network test case of study as a three-phase impedance matrix per branch for 
calculating voltage drop across the system. Finally, the seventh input is the mismatch limit that defines the 
stopping criterion, and we considered for this test a value of 1x10. 

The process section concerns all calculations required to obtain a solution of the power flow and, 
thus, the power losses of the system. We have divided the process into six steps: topology of the system, 
current injection calculation, branch current calculation (backward step), node voltage updating (forward 
step), maximum voltage deviation calculation, convergence verification, and calculation of power losses. 
There are some temporary data saved for further calculations including: the current node injection matrix; 
branch current matrix; updated voltage matrix (i.e., the same initial voltage matrix but updated with the 
branch currents calculated); and maximum voltage deviation, which is the higher difference between voltages 
in the last and current iterations. To obtain convergence, this value must be lower or equal to the mismatch 
limit. 

We divided the output section into three groups. The first group corresponds to the branch currents, 
the second group is the node voltage magnitudes and angles, and the final group is the real power loss. In the 
following subsections, we explained how BF-PF works for solving the power flow problem in radial 
distribution networks. As an example, the IEEE 13-node test feeder has been selected to develop the 
complete process. 
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Figure 1. Flowchart for the standard BF-PF 


2.1. Line segment matrix 

First, the algorithm requires the topology of the system. Thus, a line segment matrix is needed to 
verify the node connections. Table 1 presents the IEEE 13-node test feeder line matrix that can be accessed 
in [29], [30]. In that matrix, columns 1 and 2 are referred to as sending and receiving nodes of a branch, 
respectively. Column 3 introduces the length of branches. Finally, column 4 specifies the impedance 
configuration code, which provides information about the impedance matrix of each branch. Those 
configurations can be found in [29] and [30]. In addition, the number of rows must be exactly the number of 
branches that are in the system. The improvements proposed in this present paper are precisely in obtaining 
the current topology of the system considering a given line segment matrix. Those improvements are 
documented in Section 3. 
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Table 1. IEEE 13-Node Test Feeder. Line Segment Data 


From To Length (ft.) Configuration 
632 645 500 603 
632 633 500 602 
633 634 0 XFM-1 
645 646 300 603 
650 632 2000 601 
684 652 800 607 
632 671 2000 601 
671 684 300 604 
671 680 1000 601 
671 692 0 Switch 
684 611 300 605 
692 675 500 606 


2.2. Initial voltage matrix 

Initial voltages are presented in a matrix that considers the number of phases of each node. Initially, 
this matrix will contain nominal line-neutral voltages represented in real or per unit values before the 
iteration process, and updating data at each iteration. The matrix has N rows, which represent the total 
number of nodes, and four columns, which represent the line-neutral voltage per phase at each node. The 
node numeration is in column 1 and line-neutral voltages for phases a, b, and c are in columns 2, 3, and 4, 
respectively. Whether a phase is not connected the initials NA (Non-Aplicable) appears. Table 2 presents the 
voltage matrix for the IEEE 13-node test feeder. Whether a regulator is found, the downstream node would 
have the same voltage that was set in the regulator. 


Table 2. IEEE 13-Node Test Feeder. Initial Voltage Matrix 


Node Phase-A (kV) Phase-B (kV) Phase-C (kV) 
650 2,402 + 0i -1,201 — 2,08: -1,201 + 2,081 
632 2,45 — 0,1070i -1,316 —2,1291 -1,141 + 2,161i 
633 2,402 + 0i -1,201 — 2,081 -1,201 + 2,081 
634 2,402 + 0i -1,201 — 2,08: -1,201 + 2,081 
645 NA -1,201 — 2,081 -1,201 + 2,081 
646 NA -1,201 — 2,081 -1,201 + 2,081 
671 2,402 + 0i -1,201 — 2,08: -1,201 + 2,081 
692 2,402 + 0i -1,201 — 2,081 -1,201 + 2,081 
675 2,402 + 01 -1,201 — 2,08: -1,201 + 2,081 
684 2,402 + 0i NA -1,201 + 2,081 
611 NA NA -1,2+ 2,081 
652 2,402 + 0i NA NA 
680 2,402 + 0i -1,201 — 2,08: -1,201 + 2,081 


2.3. Load matrix 

A load matrix represents the power loads installed at each node per phase and it is required for 
further calculations. This matrix must have seven columns and N rows. The first column has the node 
numeration while columns 2—7 contain information concerning the active and reactive load on each node. 
They are organized in pairs per phase; for example, active and reactive load for phase a are in columns 
2 and 3, respectively, and so on for phases b and c. Table 3 can be checked in [30]. Negative numbers 
represent a reactive power injection from capacitors installed in the node. 


Table 3. IEEE 13-Node Test Feeder. Load Matrix 
Node Phase-A (kW) Phase-A (kKVAR) Phase-B (kW) Phase-B (kVAR) Phase-C (kW) Phase-C (k VAR) 
650 0 0 0 0 0 0 


632 8,5 5 33 19 58,5 34 
633 0 0 0 0 0 0 
634 160 110 120 90 120 90 
645 0 0 170 125 0 0 
646 0 0 230 132 0 0 
671 393,5 225 418 239 443,5 254 
692 0 0 0 0 170 151 
675 485 -10 68 -140 290 12 
684 0 0 0 0 0 0 
611 0 0 0 0 170 -20 
652 128 86 0 0 0 0 
680 0 0 0 0 0 0 
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2.4. Node current injection matrix 

A current injection matrix is defined as the required current by each single load installed in each 
node. For obtaining the current through the branches, calculation can be performed with (1), where n is the 
node numeration and k the kth iteration. 


JE [2 
no ve 


(1) 


For the IEEE 13-node test feeder, the matrix obtained before iterations is presented in Table 4. The 
current injection matrix exhibits N rows and four columns, where the first column contains the node 
numeration, and from the second to the fourth columns the currents per phase are reported. 


Table 4. IEEE 13-Node Test Feeder. Current Injection Matrix 


Node Phase-A Phase-B Phase-C 
Amp Amp Amp 

650 0 0 0 
632 3,374 - 2,18851 -13,391 - 7,2251 1,128 + 27,6671 
633 0 0 0 
634 66,617 - 45,7991 -57,433 - 24,5331 7,47 + 62,005i 
645 NA -80,463 - 35,2761 0 
646 NA -95,477 - 55,4531 0 
671 163,837 - 93,68 1i -173,197 - 100,9661 -0,741 + 212,793: 
692 0 0 19,057 + 92,7331 
675 201,934 + 4,164i 36,325 - 53,6641 -56,045 + 107,065i 
684 0 NA 0 
611 NA NA -42,602 + 57,134: 
652 53,294 - 35,8071 NA NA 
680 0 0 0 


2.5. Branch current matrix—backward step 
With the current injection matrix, it is possible to know the current flowing for each branch of the 


distribution network by applying (2) and (3). (2) refers to the current flowing in a branch connected from 
node 7 to a load node j. (3) refers to the current flowing in a branch connected from node h to several 
downstream nodes m that could have loads. This process is made from final nodes back to the feeder. In these 
equations, / is the current flowing from node i to node j. I; is the current injected to the load at node i. J; is 
the current injected to the load at node j. Iim is the current flowing from node i to the other nodes m. J); is the 
current flowing from an upstream node A to the node i and m are the nodes downstream connected to node i. 


Lee (2) 


Li D +I; 


(3) 


Table 5 shows the branch currents before iteration for the IEEE 13-node test feeder. The first two 
columns are used for node number identification (sending and receiving nodes). The next three columns are 
used for the current flowing through each phase. The number of rows must be equal to the number of 


branches. 


Table 5. IEEE 13-Node Test Feeder. Branch Currents 


Node from Node to Bear ee ee 
632 650 489,056 - 173,311i -383,636 - 277,117i -71,733 + 559,399 
633 632 268,551 - 41,6361 -21,109 - 78,1971 -29,518 + 261,804: 
645 632 0 -175,94 - 90,7291 0 
671 632 217,131 - 129,487i -173,197 - 100,9661 -43,343 + 269,928: 
634 633 66,617 - 45,7991 -57,433 - 24,5331 7,47 + 62,0051 
692 633 201,934 + 4,164: 36,325 - 53,664: -36,988 + 199,799 
646 645 0 -95,477 - 55,4531 0 
684 671 53,294 - 35,8071 0 -42,602 + 57,134i 
680 671 0 0 0 
675 692 201,934 + 4,164: 36,325 - 53,664: -56,045 + 107,065i 
611 684 0 0 -42,602 + 57,134i 
652 684 53,294 - 35,8071 0 0 
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2.6. Node voltage matrix—forward step 

With the forward step, it is possible to update the voltages in each node of the system. This matrix 
allows identifying the voltage profile while the branch currents are updated in the previous step as the voltage 
drop depends on the branch impedance and the current flowing on the branch. 

Node voltages are updated using both the configuration impedance matrix [29], [30] and the branch 
current matrix by using the general (4). PV nodes will have the same nominal value in all iterations, such 
node is the substation and all downstream nodes connected to a regulator. This process starts from the 
substation node toward the final nodes. 


V, =V; -Z; xl; (4) 


Equations (5)-(11) are the special conventions obtained from (4), considering that we can calculate 
the voltages for three-phase, two-phase, and single-phase branches. 


ja Va aa-n LZ ab-n Zn ij-a 
Vin = V, r Zra-n Lyn IAIN I; (5) 
[Je] L Vo J3 L^ ca-n cb-n ce-n ij-c 
Va = Va ra L ieee Z ob-n Ie] (6) 
| jb L V, J3 L^ ba-n Z pb-n Lj 
Vin 2 (V; | = | ge Lia, l ie | (7) 
Vie J [Ve 4d a= TA PR 
ARAR Z Lj 
Ja = ia | aa-n ac-n l y-a | (8) 
L je L V, J ca-n VAT Tiz 
Va = Va 7 Z wan x Eig (9) 
Vin =V p — Zob-n xI» (10) 
V, =V, -Zon Xl, (11) 


The subscripts a, b, and c refer to the branches’ phases. The impedance matrix is a reduction found 
in all IEEE test case formats [29], [30]. This reduction is derived from the Carson equations, which can be 
consulted in [31]. Table 6 presents the node voltages updated in the first iteration concerning the IEEE 
13-node test feeder. This matrix must have N rows and four columns. The first column is for node 
numeration while voltage per phase is registered in columns 2—4, considering phases a, b, and c. 


Table 6. IEEE 13-Node Test Feeder. Node Voltage Matrix 


Node Phase-A Phase-B Phase-C 

kV kV kV 
650 2,402 + Qi -1,201 — 2,081 -1,201 + 2,081 
632 2,45 — 0,107i -1,316 — 2,129i -1,141 + 2,16li 
633 2,436 — 0,134i -1,319 — 2,134i -1,117 + 2,136i 
634 2,360 — 0,172i -1,312 — 2,0681 -1,064 + 2,097i 
645 NA -1,319 — 2,108i -1,142 + 2,169i 
646 NA -1,321 — 2,102i -1,142 + 2,171i 
671 2,408 — 0,144i -1,329 — 2,091 -1,068 + 2,147i 
692 2,433 —0,222i -1,325 — 2,186i -1,055 + 2,089i 
675 2,418 — 0,232i -1,334 — 2,1871 -1,053 + 2,085i 
684 2,406 — 0,145i NA -1,064 + 2,146i 
611 NA NA -1,056 + 2,145i 
652 2,392 — 0,142i NA NA 
680 2,408 — 0,144i -1,329 — 2,090i -1,068 + 2,147i 
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2.7. Convergence requirement 

The BF-PF has an iterative process and after each single iteration the node current injection, branch 
currents, and node voltages are updated until the convergence rate is fulfilled. The minimum allowed 
mismatch in the absolute subtraction between voltage magnitude of last and current iteration is 1 x 10“ for 
each node. The error can be calculated at each iteration as shown in (12), where e is the mismatch and < is the 
current iteration. 


labs (1) — abs (V;") <e (12) 


2.8. Real power loss calculation 

When convergence is reached, the total real power losses can be computed using (13), where / is the 
branch number, rea/(Z;) is the real part of the impedance matrix of the /th branch, and /; is the current flowing 
per branch / in each phase. 


nbr 


Pg = > ,veal(Z,) xI; (13) 
[=] 


One of the main problems of programming this method is that the BF-PF requires that the nodes are 
organized in an ascendant manner. That means that a node number 2 can be fed by a node number 1, 
otherwise the algorithm will show some errors. In the next section, we present a modified part of the 
algorithm to deal with different node number orders. 


3. MODIFIED BF-PF 

In this section, we present a procedure applied to identify the topology and calculate branch currents 
and node voltages in a radial unbalanced distribution network, independent from the node number order. 
As in the standard BF-PF, the inputs considered in the algorithm are given by the feeder test case. The 
outputs are the same as those presented in Figure 1, namely, branch currents, nodes voltages, and active 
power losses. 


3.1. Loading the case study 

The first step is to load the feeder test case that contains the slack node, the number of switches, and 
the data of lines, nodes, and regulators. The number of branches is expressed in (14) and can be used to 
verify that in case of multiple changes in the topology, the new network is radial, where nbr is the number of 
branches and N is the total number of nodes. 


nbr = N-1 (14) 


3.2. First radiality constraint verification 

In this step, we verify that the number of nodes is the same to that calculated in (14) when the data 
from the feeder test case is obtained. For example, in the IEEE 13-node test feeder, the number of nodes is 
N=13 and the number of branches is nbr=12. In this example, the algorithm must identify if a low number of 
branches is connected and, thus, the distribution network does not comply with radiality. 


3.3. Neighborhood matrix 

The neighborhood matrix is defined as the identification of neighboring nodes [4]. This approach 
shows a list of connections between nodes with the rest of the system. As an example, the neighborhood 
matrix for the IEEE 13-node test feeder is presented in Table 7. It must contain two columns and two times 
nbr rows. Column 1 presents the node number and column 2 refers to all nodes connected to the node in 
column 1. If a node is connected to more than one other node, then it appears more than once in the first 
column. From Table 7, node 632 is connected to 650, 633, 645, and 671. 


3.4. Determining the type of node 

Once the neighborhood matrix is built, a mf mt matrix is necessary to identify the number of 
branches connected to each node. Table 8 presents the mf_mt matrix for the IEEE 13-node test feeder, which 
has N rows and four columns. The node number is in the first column. The starting row, where node n is 
found into the first column of the neighborhood matrix, is placed in column 2. The ending row, where the 
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node n is in the first column of the neighborhood matrix, is specified in column 3. Column 4 contains the 
subtraction between columns 3 and 2 [27]. 


Table 7. IEEE 13-Node Test Feeder. Neighboring Matrix 
Column 1l-part 1 Column 2-part 1 Column |—part 2 Column 2-part 2 


Node Connected to Node Connected to 
650 632 671 680 
632 650 671 684 
632 633 671 692 
632 645 680 671 
632 671 684 671 
633 632 684 652 
633 634 684 611 
634 633 652 684 
645 632 611 684 
645 646 692 671 
646 645 692 675 
671 632 675 692 


Table 8. IEEE 13-Node Test Feeder. Mf_mt Matrix 


Node number mf (initial row) mt (ending row) mt-mf- 
650 1 1 0 
632 2 5 3 
633 6 7 1 
634 8 8 0 
645 9 10 1 
646 11 11 0 
671 12 15 3 
680 16 16 0 
684 17 19 2 
652 20 20 0 
611 21 21 0 
692 22 23 1 
675 24 24 0 


As shown in this table, the node 650 starts in the first row of the neighborhood matrix and ends in 
this same row, which indicates that it has only one connection through a branch to one other node (632). 
Moreover, node 632 starts in row 2 from the neighborhood matrix and ends in row 5, which indicates that 
this node is connected through four branches to four nodes. 

Column 4 presents the case when a node is an ending node, intermediate node, or a junction node. 
This is determined as follows [27]: 

1. An ending node refers to the node that connects only to another node through one or more branches. For 
this case, the value in column 4 is zero. 

2. An intermediate node can be identified as the node connected to two other nodes through one or more 
branches for each pair of nodes. For this case, the value in column 4 is 1. 

3. A junction node is the element that is used to connect more than two nodes with more than two branches. 
For this case, the value in column 4 is greater than 1. 

The slack node must be found in the mf_mt matrix and identify its position, and not to be considered 
as an ending node. Because the slack node could have one or more connections, it must be considered as the 
root when searching the topology of the system. The binary search initializes from the slack node to find each 
connection to the rest of the nodes; therefore, we considered in this paper to explain how we deal with the 
topology to avoid mistakes when the connections change, and the search is needed to reconfigure. 


3.5. Topology testing 

A verification process is done to identify the network topology and constraints such as radiality. 
First, two matrices are defined to identify the node connections. Those are Line _ prove and Node _ prove 
matrices. The Line prove matrix registers, one by one, the connection between the nodes. The Node_prove 
matrix identifies when a node has been included in the Line prove matrix, which means that we have 
encountered the upstream connection of that node. 

The Line_prove matrix contains nbr rows and two columns. Column | contains the receiving node 
and column 2 the sending node per branch. Although, this matrix is initiated full of zeros, the data must be 


Int J Elec & Comp Eng, Vol. 9, No. 1, February 2019 : 85 - 101 


Int J Elec & Comp Eng ISSN: 2088-8708 Oo 93 


updated iteratively. The Node_prove matrix has N rows and two columns. Column | is referred to as the node 
numeration while column 2 is a binary value assigned to each node that allows identifying if that node was or 
not evaluated before and is initially full of zeros. 

The first step is finding the slack node and its branches. If the slack node is connected to only one 
node, then the first row of the Line_prove column will contain in column | the node that is downstream 
connected to the slack node and the slack node will be in column 2. In the case where the slack node has 
more than one connection, there must be as many rows with the slack node number as there are connections. 

Once this first step is done, the node(s) connected to the slack node must be found in the 
Node_prove matrix to place the number 1 in column 2 which is interpreted as we have evaluated a node and 
found its feeding previous node. As a brief example, consider the IEEE 13-node test feeder and both its initial 
Line_prove and Node_prove matrices in Tables 9 and 10, respectively. From Table 9 and Table 10, it is 
specified that the node 632 connected to the slack node (650) is checked in the Node_prove matrix, putting 
the number | next to its position in the matrix. In this instance, only one branch has been identified. 


Table 9. IEEE 13-Node Test Feeder. Table 10. IEEE 13-Node Test Feeder. 
Initial Line_ Prove Matrix Initial Node_Prove Matrix 
Receiving node Sending node Receiving node Binary search 
632 650 632 1 

645 0 
675 0 


3.6. Filling Node_prove and Line_prove matrices 

As a Starting point, a parameter called num lines has to be created with the value of 1. This 
parameter allows checking the branches that follow the slack connections. An iterative process is initiated 
from | to nbr according to the index i that locates the last row that has been filled in the Line_prove matrix. 
From the ith row of Line_prove matrix, the receiving node is identified to find all its connections from the 
neighborhood matrix. Those connections are saved in the first column of a new variable called post_nodes. A 
second column is required to recognize which of those connections is the upstream node. This is possible 
comparing the post_nodes matrix against the second column of the Line_prove matrix, when a common node 
is encountered, the number | is placed in the second column of post_nodes next to the node identified. This 
process can be observed in Table 11. 


Table 11. IEEE 13-Node Test Feeder. Example of Post_Nodes Matrix for Node 632 


Connections Binary decision 
650 1 
645 0 
633 0 
671 0 


We can determine that the node 650 is connected to the node 632 because that is the only common 
node between the first column of post_nodes matrix and the second column of Line_prove matrix. We can 
also determine that the nodes 645, 633, and 671 are the downstream nodes connected to the node 632. Those 
connections are reported in Line_prove matrix and the number | must be updated next to those nodes in the 
Node_prove matrix. The parameter num_lines is updated by adding the amount of branches included in the 
Line_prove matrix. This process ends when num_lines is equal to nbr. 


3.7. Second radiality constraint verification 

To accomplish the second radiality constraint, we identify the network topology by using the 
Line_prove matrix and creating an adjacency matrix (A) that through a binary representation, presents the 
nodes that are connected to a branch. This new matrix has nbr rows and N columns and the initial values are 
defined as zeros. Each row in the matrix represents each registered branch from the Line_prove matrix and 
the columns contain all node numbers. A branch is always connected to two nodes; that is, only two cells 
with numbers different than zero per row. The sending node will have the number 1 and the receiving node 
the number —1. To calculate the determinant of that matrix, it must be rectangular, and it has one additional 
column if compared against the number of rows, so the slack node column must be eliminated. Once this 
matrix is built, the second constraint defined as the network radiality is verified using (15). 
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det (A) =+1 (15) 


3.8. Load matrix 

This matrix contains the loads installed at each node as explained in Section 2. This matrix contains 
four rows and N columns. The node numeration can be found in the first row and the installed load per phase 
is in rows 2, 3, and 4, respectively. 


3.9. Impedance branch matrix 

The impedance branch matrix is used for updating voltage matrix and calculating real power losses 
of the radial system. Based on the configuration code, there are different 3 by 3 matrices for the three-phase, 
two-phase, and single-phase branches. From line segment data, each branch has a length value and a 
configuration code. With this information, an impedance matrix is created with three columns and three times 
nbr rows. Each branch will have three rows and three columns, constituting a 3X3 independent matrix. 


3.10. Parameter setting 

The maximum difference between voltage magnitudes from the last iteration to the current one is 
initially stated as a big number (higher than e). For this work, it has been considered as 99. Then, the 
mismatch limit is set as e=1x10~ and the number of iterations is set to zero. While the convergence criterion 
is not accomplished, the BF-PF updates the node injection currents, branch currents, and node voltages, 
during each iteration. 


3.11. Node current injection matrix 

As explained in Section 2, it is the required current per node as calculated using (1). This matrix 
contains four rows and N columns. The node numeration is in the first row and the required current per phase 
is in rows 2, 3, and 4, respectively. Each column contains the information about each node. 


3.12. Backward step-branch current calculation 

Although this step was explained in Section 2, we modified the algorithm for calculating branch 
currents, regardless of the number of nodes. To make such a calculation, it is necessary to use the Line_prove 
matrix and create a branch current matrix that contains only one column, and three times nbr rows. Each 
branch has three rows and one column, considering the current per phase. The initial values of the matrix are 
zeros and the following process is based on (2) and (3). 

There are two iterative processes, with one cycle dependent on the other: this is known as “for cycle 
inside another for cycle” in programming. The first cycle is followed with the aid of an index i, which is set 
from 1 to nbr. The order of branches is the same as they are in the Line_prove matrix. A variable called x is 
initialized at zero, which guarantees that the backward step is executed from the last node to the slack node, 
namely, from the end to the beginning, as a backward sweep. This is achieved because the x variable value is 
subtracted from nbr; as the x value is increased per iteration, the difference between nbr and x is reduced, 
thereby attaining an evaluation of each branch for the last row of the Line prove matrix to the first row. 
The branch current for the branch nbr-—x is first calculated as the reported current injected in the node current 
injection matrix in the receiving node connected to that branch. 

If node n is downstream connected to more than one node, then the second cycle is needed. 
This cycle is inside the first cycle and is accomplished with the aid of an index j, which is set from 1 to the 
number of times this receiving node is registered in the second column of the Line_prove matrix. The nodes 
in the first column of the Line_prove matrix, next to the node n, are referred as its downstream connected 
nodes and the current injected to those nodes is added to the branch current calculated in the first cycle. 
When the above process is performed, the variable x is updated as the value of the index i, starting all over 
again until į is equal to nbr. Finally, each magnitude of a branch current is calculated from the branch current 
matrix that is in complex numbers. The reason to obtain the magnitude of a branch current is for calculating 
real power losses. 


3.13. Forward step-node voltage matrix calculation 

At this step, two assumptions are considered in the algorithm. The first assumption is that the slack 
node and voltage regulated nodes keep the same voltage at every iteration. The second assumption is that the 
Line_prove matrix is used to follow a path from the slack node to the final nodes. The line segment matrix is 
also required for checking the configuration per branch, which is important data to calculate voltage drops 
between nodes using (5)-(11). An index i leads a cycle comparing the information saved in the Line_prove 
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matrix against the line segment matrix to obtain the number of phases of each branch. Depending on the 
connected phases, the algorithm must use one of (5)-(11). This algorithm saves two voltage matrices: the last 
iteration voltage matrix and the current iteration voltage matrix. 


3.14. Convergence requirement 

The convergence requirement is defined because a minimum error must be reached to determine if 
the calculations are correct and gives a stop criterion to the algorithm to finish all the process. From last step, 
two voltage matrices were obtained from last iteration and current iteration. The absolute value of the 
magnitude difference between each pair of rows of those two matrices is calculated and saved in a matrix 
called the “mismatch matrix.” The maximum encountered value in this matrix is compared with the 
mismatch limit e as presented in (12) and set as 1x10. If convergence is not reached, then the voltage matrix 
for the last iteration is updated with the values of the voltage matrix at the current iteration and the resulting 
matrix is finally ready to be used in the next iteration; otherwise, the algorithm finds the final convergence 
and the algorithm stops. 


4. TEST FEEDERS FOR RECONFIGURATION IMPLEMENTATION 

There are four distribution test cases developed by Kersting for the IEEE in [29] and [30] 
(IEEE 13, 34, 37, and 123-node distribution test feeders). In this work, we have selected the smallest test 
feeder (13-node test feeder) and largest test feeder (123-node test feeder) and modified them by adding 
switches connected to some nodes. Figure 2 and Figure 3 show the modifications applied to the two test cases 
based on [32] and [33]. 


—— 650 


646 645 632 633 634 
e 123 + 123 12s * 123 35 + 


652 680 
m Tie switch 
—— Sectionalizing switch 
Red numbers Phases 
È Regulated node 


Figure 2. Modified IEEE 13-node test feeder 


From Figure 2 and Figure 3, it can be observed that there are tie and sectionalizing switches. A tie 
switch is a normally opened switch, and all must be three-phase to allow for a reconfiguration process; 
meanwhile, a sectionalizing switch is a normally closed switch. In both figures, the numbers in red refer to 
the phases per branch, where 123 means that phases a, b, and c are connected through the branch, where 12 
means that phases a and b are connected through the branch, and so on. The regulated nodes have a regulator 
installed that is capable of maintaining constant voltage in that node. 

Information about the tie switches, added to the distribution networks for the 13-node test feeder and 
the 123-node test feeder, are presented in Table 12 and Table 13, respectively. Lengths and configuration of 
switches were calculated using the distance to the neighboring nodes. Sectionalizing switch lengths and 
configurations are already in the line segment matrix of each case study, although the nodes where they are 
connected are in Table 14 and Table 15 for the IEEE 13-node test feeder and IEEE 123-node test feeder, 
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respectively. For the case of the IEEE 123-node test feeder, the first columns in Tables 13 and Table 15 
indicate the number of the switch. 


—_— Tie switch 
Ss Sectionalizing switch 


Rednumbers Phases 


D Regulated node 


Figure 3. Modified IEEE 123-node test feeder 


Table 12. IEEE 13-Node Test Feeder. Table 14. IEEE 13-Node Test Feeder. 
Tie Switches Line Segment Information Sectionalizing Switches 
Fro To Length Configuration From node To node 
m node (f) code 632 671 
node 671 692 
646 684 2000 601 
633 692 2000 601 
Table 13. IEEE 123-Node Test Feeder. Table 15. IEEE 123-Node Test Feeder. 
Tie Switches Line Segment Information Sectionalizing Switches 
Switch From To Length Configuration Switch number From node To node 
number node node (ft) code 1 18 21 
16 56 92 350 6 2 23 25 
17 83 95 1975 2 3 25 28 
18 25 35 906 2 4 42 44 
19 250 300 1364 3 5 57 60 
20 51 65 600 4 6 67 72 
21 101 151 1600 3 7 80 81 
22 79 450 1304 5 8 86 87 
9 101 105 
10 13 152 
11 18 135 
12 60 160 
13 97 197 
14 54 94 
15 151 300 
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5. RESULTS AND DISCUSSION 

In this section, we present the results for different combinations of switches in both the IEEE 
13-node and IEEE 123-node test feeders to demonstrate the effectiveness of the proposed modifications to 
the BF-PF. For the smallest case, all possible combinations of switches are recorded. However, for the largest 
case, the number of possible combinations is 4,194,304 and feasible combinations number 497,420; thus, we 
listed only the best 20 combinations. 


5.1. Results for the modified IEEE 13-node test feeder 

The modified IEEE 13-node test feeder has four switches, indicating that the possible number of 
combinations is 16 as the expression for obtaining the number of combinations in a binary decision variable 
is calculated from (16), where n indicates the number of switches installed in the system [34]. 


comb = 2” (16) 


Table 16 shows the results for all possible combinations. The opened switches for every case are in 
column 1 while the real power losses are in column 2. If a combination is not feasible, then the algorithm 
displays a message that allows us to know exactly what constraint is not fulfilled. 


Table 16. IEEE 13-node test feeder. Evaluation of the Modified BF-PF 


Opened switches Losses (kW) 

All Constraint: Minimum number of branches 
646-684, 633-692, 671-692 Constraint: Minimum number of branches 
646-684, 633-692, 632-671 Constraint: Minimum number of branches 
646-684, 633-692 102,98 

646-684, 671-692, 632-671 Constraint: Minimum number of branches 
646-684, 671-692 73,61 

646-684, 632-671 125,78 

646-684 Constraint: Minimum number of branches 
633-692, 671-692, 632-671 Constraint: Minimum number of branches 
633-692, 671-692 Constraint: Not feeding all nodes 
633-692, 632-671 109,83 

633-692 Constraint: Minimum number of branches 
671-692, 632-671 68,43 

671-692 Constraint: Minimum number of branches 
632-671 Constraint: Minimum number of branches 
None Constraint: Minimum number of branches 


Five combinations could fulfill radiality constraints and reach the convergence requirement. 
The best combination for minimal real power losses is when the switches between nodes 671—692 and nodes 
632-671 are opened, obtaining a total real power loss of 68.43 kW. This test feeder is easy to evaluate 
because it does not require large computational resources to calculate the power flow for the small number of 
combinations. 

The unfeasible combinations are separated into two groups. The first group corresponds to the 
configurations that obtain the number of branches different than calculated in (14), preventing to perform the 
power flow. The second group corresponds to the configurations that do not feed all nodes, which is the case 
when switches between nodes 633—692 and between 671—692 are opened, and node 675 is isolated. For this 
case study, the modified BF-PF can accomplish good results, allowing its use for small-scale radial 
unbalanced distribution networks. 

Aditionally, the five feasible combinations obtained with the IEEE 13-node test feeder were 
compared to validate the effectiveness of the Modified BF-PF (MBF-PF) against a Standard BF-PF 
(SBF-PF) [27]. Such comparison is presented in Table 17. 


Table 17. Comparison between MBF-PF and SBF-PF (IEEE 13-Node Test Feeder) 
Opened switches MBEFE-PF response SBF-PF response 


646-684, 633-692 Converge Does not converge 
646-684, 671-692 Converge Does not converge 
646-684, 632-671 Converge Does not converge 
633-692, 632-671 Converge Does not converge 
671-692, 632-671 Converge Does not converge 
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From Table 17, it is clear that the SBF-PF does not find a solution for network reconfiguration, 
because it needs the special numeration considered in the MBF-PF. This result probes that the MBF-PF can 
find solutions for reconfiguration regardless the numeration of nodes. 


5.2. Results for the modified IEEE 123-node test feeder 

The modified IEEE 123-node test feeder is a large-scale radial unbalanced distribution network. 
From Table 13 and Table 15 it is observed that there are 22 switches. Hence, 13 switches must be closed and 
9 switches must be opened to fulfill radiality constraints. There are 4,194,304 combinations to evaluate, but 
to obtain only the feasible combinations, which number is 497,420. (17) can be used [34], where perm is the 
number of permutations, n the total number of switches, and r the number of whichever required opened or 
closed switches (as this equation is symmetric). The symbol (!) refers to the factorial value. 


n! 


pst r\(n—r)! 


(17) 


Despite the feasible combinations is 11.86% of the total possible combinations it is still a very large 
effort to evaluate all combinations in reasonable time for a real operation system case. Reconfiguration is a 
challenge to large-scale radial unbalanced distribution networks. 

We evaluated all feasible combinations using the proposed algorithm in an AMD A10-5745M APU 
with Radeon(tm) HD Graphics 2.10 GHz, 8 GB RAM computer in 36,022.39 seconds (approximately 10 
hours). The best 20 combinations are reported in Table 18, which contains the 22 switches listed as they are 
in Table 13 and Table 15. When a switch is closed, the number | appears below the number of that switch, 
otherwise 0. The real power losses for each configuration are shown in the last column. 


Table 18. IEEE 123-Node Test Feeder. Evaluation of the Modified BF-PF. Best 20 Combinations 
State of switches: 0 for opened, | for closed Losses (kW) 
9 10 11 12 #13 14 15 16 
204,16 
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From Table 18, it is observed that the best configuration for minimizing the real power losses is to 
keep the switches from nodes 67 and 72, 80 and 81, 151 and 300, 56 and 92, 25 and 35, 250 and 300, 51 and 
65, 101 and 151, 79 and 450 opened. When an unfeasible combination is evaluated, the algorithm prompts 
the same message errors as those reported in Table 16. The worst scenario presented in Table 18 is 2.28% 
higher than the best combination obtained with the complete evaluation performed for all feasible 
combinations. This result means that the algorithm is capable to find good solutions with least amount of 
combinations leaded by a random process bounded by natural processes. Aditionally, the comparison 
between MBF-PF and SBF-PF [27] shows that the SBF-PF could not find the solutions to the same 
combinations, because it requires to renumber all the nodes for each evaluation proposed for the network 
reconfiguration. 


6. CONCLUSION 
The standard BF-PF is a known algorithm for evaluating power flow in distribution systems. Some 
tools developed by researchers follow a strict order in node numeration. If that numeration order is changed, 
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then the standard BF-PF is not capable to reach the convergence requirement or even evaluate the power 
flow. To avoid this complication, the authors have proposed a modified BF-PF algorithm that carries a binary 
search from the substation node to the last node in the system. This is accomplished with the aid of several 
matrices that guarantee that the substation or slack node will always be the root of the tree search process, 
obtaining a path whereby the current is flowing. 

This purpose was proven on two known benchmarks from the IEEE: the IEEE 13-node test feeder 
and the IEEE 123-node test feeder. With these evaluations, the authors could demonstrate that with these 
developed modifications, this algorithm becomes a powerful tool to tackle both small- and large-scale radial 
unbalanced distribution networks, regardless the numeration of the nodes. 

Modifications on case studies were presented to evaluate different topologies in each test feeder, 
considering both feasible and unfeasible combinations. The algorithm was capable to evaluate every single 
combination and, when an unfeasible combination was evaluated, a message prompted showing what 
constraint was not fulfilled by the configuration. The authors also introduced a permutation equation for 
obtaining only feasible combinations. This allows making power flows in a reconfiguration process while 
avoiding unfeasible combinations. The results for each case study were presented and analyzed to comment 
the effectiveness of the proposed modifications. 
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